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We study dynamical phases of a driven Bose-Einstein condensate coupled to the light held of a 
high-Q optical cavity. For high held seeking atoms at red detuning the system is known to show a 
transition from a spatially homogeneous steady-state to a self-organized regular lattice exhibiting 
super-radiant scattering into the cavity. For blue atom pump detuning the particles are repelled 
from the maxima of the light-induced optical potential suppressing scattering. We show that this 
generates a new dynamical instability of the self-ordered phase, leading to the appearance of self- 
ordered stable limit-cycles characterized by large amplitude self-sustained oscillations of both the 
condensate density and cavity held. The limit-cycles evolve into chaotic behavior by period doubling. 
Large amplitude oscillations of the condensate are accompanied by phase-slippage through soliton 
nucleation at a rate which increases by orders of magnitude in the chaotic regime. Different from a 
superhuid in a closed setup, this driven dissipative superhuid is not destroyed by the proliferation 
of solitons since kinetic energy is removed through cavity losses. 



FIG. 1. A Bose-Einstein condensate of neutral atoms (blue 
surface) trapped inside an optical resonator is driven with a 
monochromatic laser with amplitude Q. The laser frequency 
is blue-detuned by an amount A a with respect to an internal 
atomic transition g o e and by A c with respect a standing- 
wave cavity mode ~ cos (k c x) dehning the characteristic recoil 
frequency cjr = hk c /2m with atomic mass m. The coupling 
between a single atom and the cavity mode has a strength 
go. For large A a (A c , ft, go, cjr), with k being the cav¬ 
ity linewidth due to leakage out of the mirrors, the dispersive 
regime is reached, where the atoms stay effectively in their in¬ 
ternal ground state. They experience negligible spontaneous 
emission and an optical potential (orange surface)arising from 
the interfering cavity and pump fields (see Eq.Q). For sim¬ 
plicity we assume a sufficiently strong extra atom trapping 
in the radial y , z directions to confine the dynamics onto the 
cavity axis x and generate an effective ID geometry. 


I. INTRODUCTION 

The experimental realization of strong collective cou¬ 
pling between a large number of (ultra)cold atoms and 
the electromagnetic field of Fabry-Perot cavities w, 


nanophotonic fibers pm or photonic crystals [12] opens 
up new interesting routes both in the field of quantum op¬ 
tics and condensed matter physics. From the latter point 
of view, the interesting new ingredient is provided by 
very well controllable strong long-range photon-mediated 
atom-atom interactions appearing due to the back-action 
of even a single atom onto the light field. Strong long- 
range interactions can indeed lead to several intriguing 
phenomena mm and are crucial in many intensively 
explored condensed matter phases, like supersolids m 
or topological states m- What is more, differently from 
typical condensed matter situations, these light-mediated 
interactions are in general i) retarded, since the photon 
field owns intrinsic timescales which can be made compa¬ 
rable with atomic scales, ii) non-conservative, since the 
system is typically driven and dissipates energy through 
atomic spontaneous emission and photon losses through 
the cavity mirrors. 

A striking consequence of cavity field mediated in¬ 
teractions is the appearance of self-ordered phases [ 13 - 
124] where the particles break a translation symmetry by 
forming a spatial pattern determined by the characteris¬ 
tic interaction length scale. This phenomenon has been 
observed experimentally both for a thermal gas mm and 
an ultracold Bose-Einstein condensate (BEC) [5j \7\ cou¬ 
pled to a standing-wave mode of an optical cavity as 
sketched in Fig. [l] In the regimes considered so far, a 
thermal gas and a BEC share the same qualitativ be¬ 
havior |25]. Interestingly the self-ordering of a BEC can 
be closely mapped to the superradiant transition of the 
famous Dicke model of N two-state atoms coupled to a 
single cavity mode Ham]. If fundamental collective me¬ 
chanical excitations of the BEC are effectively treated as 
quantum mechanical oscillators, one also gets a simulator 
for a general optomechanics setup close to zero temper¬ 
ature [23128] . 

Here we study a new regime of quantum gas cavity 
QED [29], where the self-ordering phase transition is 
tight to dynamical instabilities. This gives rise to new 
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types of nonequilibrium phases for which the peculiarities 
of the BEC play a crucial role. This novel behaviour ap¬ 
pears upon a rather innocent looking change of operating 
conditions, namely by choosing the frequency of the driv¬ 
ing laser larger than the atomic internal transition fre¬ 
quency (blue detuning). In this regime the atoms are low 
field-seekers as opposed to high field seekers in the typi¬ 
cally considered red detuning case (see Fig. [l]). Naively, 
one would expect that this prevents any self-ordering as 
the atoms are pushed towards field minima, where light 
scattering is suppressed. Surprisingly, a closer look re¬ 
veals that the complex interplay of collective coherent 
scattering and optical dipole forces still can generate a 
self-ordered phase at sufficient pump strength. How¬ 
ever, the particles are now localized at cavity field nodes 
and this order gets dynamically unstable again at only 
a somewhat higher critical pump intensity, as illustrated 
in Fig. [2] 

Interestingly, this instability does not simply lead to 
heating and disintegration of the order, but we find the 
emergence of limit-cycles, whereby the condensate per¬ 
forms large periodic self-sustained oscillations between 
different ordered patterns. The atomic density oscilla¬ 
tions are tightly coupled to the oscillation of the cav¬ 
ity field with the same frequency, as shown in Fig. Bv 
This provides a build-in non-destructive monitoring tool 
of the nonlinear dynamics. By further increase of the 
drive strength the limit-cycles turn into chaotic dynam¬ 
ics by doubling their period (see Fig. [3|. 

Dynamical instabilities toward limit-cycles evolving 
into chaos have been also observed with nanomechani¬ 
cal oscillators coupled to light [301432] . Limit-cycles have 
been studied within the open Dicke model as well 33], 
where chaos appears in the closed-system limit [34h36] . 
Self-sustained oscillations of a BEC inside a driven cavity 
have been observed and shown to be very well de¬ 
scribed through an optomechanical model. In the same 
setup, a transition to chaos has also been predicted [38] . 
analogous to the one appearing with nonlinear dielectric 
medium [39]. Recently, superfluid Josephson dynamics 
of a BEC coupled to a driven cavity have also been the¬ 
oretically studied m • 

In this work, we show how the dynamical self-ordered 
regimes of large nonlinear excitations let the BEC hall¬ 
marks emerge clearly, which prevents an understanding 
of the system via (generalized) optomechanical models. 
In particular, once direct short-range atom-atom inter¬ 
actions are taken into account, the superfluid nature of 
the BEC manifests itself through the onset of phase- 
slippage dynamics mug, whereby the condensate low¬ 
ers the kinetic energy stored in the phase of the macro¬ 
scopic wave function by creating phase singularities in 
the form of nonlinear dispersive waves (solitons in our 
one-dimensional model). While phase-slips take place 
periodically in the lymit-cycle phase, they appear irreg¬ 
ularly and at a much faster rate in the chaotic regime, 
as illustrated in Fig. [4] Interestingly, different from a su¬ 
perfluid in a closed system where phase-slip proliferation 


eventually destroys the superfluid, here the dissipation 
through cavity losses counteracts this heating process by 
subtracting energy from the system. 

Our findings introduce a new scenario where nonlinear 
chaotic dynamics, self-organisation and superfluidity ap¬ 
pear together in a driven/dissipative system, bridging be¬ 
tween the optomechanics/nonlinear optics and condensed 
matter communities. 


II. DYNAMICAL PHASE DIAGRAM 



FIG. 2. A) Nonequilibrium phase diagram for the driven 
Bose-Einstein condensate as in the setup of Fig. [l] as a func¬ 
tion of the effective pump strength g = y/Ngo^l/ A a and the 
detuning A c , for k, = IOcjr, UqN = 12.1o;r, and g aa = 0. The 
phase is determined by the long-time behavior of the system 
starting from an empty cavity and a spatially homogenous 
condensate along x. The color scale indicates the growth rate 
Rea; of the most unstable collective excitation mode above 
the steady-state. The corresponding time evolution of the 
cavity mode amplitude a is shown in panels I,II, and III. 
Four different phases emerge. In the normal (N) phase, the 
steady-state is still an empty cavity with a homogeneous con¬ 
densate. In the self-organised (S-O) phase, the steady state is 
stable and has a finite a accompanied by the corresponding 
condensate density modulation in space, optimising the light 
scattering into the cavity. In the self-organised limit-cycle (S- 
O-L-C) phase, the steady state is unstable and evolves into 
periodic self-sustained oscillations of large amplitude about 
a finite value of a. At every time during the oscillation the 
condensate is self-organised, i.e. has chosen the density mod¬ 
ulation giving rise to the instantaneous cavity amplitude. In 
the self-organised chaotic (S-O-Chaos) phase, the collective 
oscillations loose their periodicity. The transition from limit- 
cycles to chaos takes place by period doubling, as illustrated 
in Fig. [3] 

As sketched in Fig. [TJ we consider an ultracold atomic 
gas trapped along the axis a single mode of a Fabry-Perot 
cavity. The gas is illuminated with coherent light which is 
blue detuned to an atomic transition. Corresponding ex¬ 
periments with BECs of ultracold atoms mm have shown 
to be well described by coupled classical field equations 
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p~8j describing the mean-field dynamics of the BEC with 
the Gross-Pitaevskii equation [45] 


ihd t 'ip(x , t) = 


2m 


+ 9aa\ip(x,t)\ + frU 0 \a(t)\ 2 cos 2 (k c x) 


+2h(r]/VN)Re[a(t)] cos (k c x) 


V>(M) (1) 


and the cavity field by its coherent component dynamics 
a(t) : 




Uo Jdx\'ip(x^ t)\ 2 cos 2 (k c x) 
{ri/y/N) Jdx\'ip(x, t)\ 2 cos(k c x) . 


— A c — IK 
+ 


a(t) 

( 2 ) 


In the Gross-Pitaevskii framework, all N particles in 
the condensate, i.e. occupying the same quantum state, 
share the same BEC wavefunction ^(x^t) normalized to 
N. The first term in Eq. 0 describes the kinetic en¬ 
ergy of the quantum motion of the particles with mass m, 
which we restricted along the cavity axis x upon assuming 
additional trapping in the other directions. This can be 
easly achieved in ultracold atom experiments [46]. Cor¬ 
respondingly, the direct atom-atom interaction strength 
g aa is the effective coupling for the one-dimensional prob¬ 
lem [45] . Uo = go/A a is the potential depth per photon 
felt by an atom (see Eq. 0 ) as well as the energy shift 
of the cavity resonance per atom (see Eq. 0). It re¬ 
sults from emission and absorption of a photon from and 
into the cavity mode and is therefore oc cos 2 (k c x) in the 
dispersive regime where the excited atomic state can be 
adiabatically eliminated. The further terms containing 
rj/y/N = go&/A a result from absorption from the pump 
and emission into the cavity mode (or viceversa) and in¬ 
troduce a further optical potential oc cos (k c x) for the 
atoms and an effective pump term for the cavity field. 
Finally, the loss of photons through the cavity mirrors 
is reflected by the field damping term — ik. Eqs. 0,0 
describe the dynamical back-action between atoms and 
light, since the optical potential felt by the atoms is de¬ 
termined by a(£), whose dynamics in turn depends on the 
atom wave function. Equivalently, the field a mediates 
infinitely long-range and retarded interactions between 
the BEC atoms. 

The phase diagram of Fig. [2] is obtained by solv¬ 
ing Eqs. 0,0 and analysing their long-time behavior. 
As initial conditions we choose a homogeneous BEC: 
^(x, 0) = \fn with n = N/L the system’s density in 
one dimension and an infinitesimally occupied cavity 
a( 0) <C 1, which is needed as a seed since Eqs. 0,0 
don’t include noise. A transition toward stable self¬ 
ordering in observed at a critical pump strength [47] 
h 2 r] 2 rit = (^R + 2g aa n)(S c + U 2 /5 C )/2, with the disper- 
sively shifted cavity detuning 5 C = — A c + UoN/2. This 
transition is due to the fact that above a critical pump 
strength the homogeneous condensate is unstable with 
respect to density modulations at the cavity wavelength 


A c . The latter indeed tend to optimise light scattering 
into the cavity which in turn enhances the density mod¬ 
ulation in a runaway process. This is stabilized by losses, 
providing convergence toward a steady state with finite 
cavity amplitude a ^ 0 and density modulation. This 
happens both for high (Uo < 0) and low (Uo > 0) field¬ 
seeking atoms. However, as we show here, in the latter 
case the stable self-ordered steady state exists only up to 
a second critical pump strength, above which it becomes 
dynamically unstable and no time independent steady 
state can be found anymore. This originates from the 
competition between the cavity oc cos 2 (k c x) generated 
potential pushing the particles to the cavity field nodes 
where scattering is suppressed and the pump-cavity in¬ 
terference component oc cos (k c x) pushing the particles to 
competing wavelength order. 


III. SELF-ORDERED LIMIT-CYCLES AND 
CHAOS 

Above the second critical driving strength the sys¬ 
tem resolves the above competion by performing self- 
sustained oscillations (limit-cycles) between self-ordered 
configurations of the condensate (see Fig. 4l Al) and cav¬ 
ity field, shown in Fig. [2 II, [3] A)-B), and|4]A2). 

As illustrated in Fig. 3] the limit-cycles involve sev¬ 
eral well defined frequency components, whose number 
is doubled before the system dynamics turns eventually 
chaotic. The onset of chaos by period doubling is a 
well known phenomenon in nonlinear systems and has 
been observed with nanomechanical oscillators coupled 
to light m- Here, like in these systems, chaos develops 
at zero temperature and without external periodic modu¬ 
lations or (delayed) feedback control. The characteristic 
timescale is given by the recoil frequency cjr which is an 
intrinsic property of the BEC-cavity system. 

Different from what so far studied with nanomechan¬ 
ical media, nonlinear dielectrics and non-directly-driven 
BECs [37, [38], the oscillations found here have a large 
(non-perturbative) amplitude where the light intensity 
scales with N due to the self-ordering of the atoms. The 
latter generates non-trivial time-dependent spatial struc¬ 
tures. 


IV. DRIVEN-DISSIPATIVE PHASE-SLIPPAGE 

While spatial structures with a length scale set by the 
cavity wavelength A c would be present had we consid¬ 
ered any (e.g. thermal) ensemble of driven polarisable 
particles, the peculiar nature of their size and dynamics 
that we oberve here originates from the combination of 
macroscopic phase coherence (encoded in the BEC wave- 
function) together with short-range atom-atom interac¬ 
tions (inducing the term gaalV’l 2 term in Eq. 0 ). These 
two properties in turn are at the core of the superfluid 
behavior of the BEC. They indeed provide the conden- 
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FIG. 3. Transition from limit-cycles to chaos. The Fourier transform of the cavity amplitude time-oscillations (left column) 
is shown together with the condensate wave function recurrence 33 2 ) = f dxj / ip(x,t 1 ) — ip(x,t 2 )\ 2 (middle column) 

and phase-space trajectories in the \a\ — |'0(fc c )| plane (right column). Here ip(k c ) is the Fourier component of the condensate 
wavefunction calculated at the cavity wave number k c . For the same parameters as in Fig. [2] results are shown for three different 
pump strengths: A) rj = 9.5cjr, B) 77 = IOcjr, and C) 77 = 11c<;r at A c = 4c<jr. The spectra clearly show period doubling 
(between A) and B)) before the onset of the noisy background in C). This indication for the appearance of chaos is confirmed 
by the large-scale structures in the recurrence R , as opposed to the stripes characterising a periodic behavior, visible in A) and 
B). In the latter, period doubling shows up as further small-scale structures perpendicular to the stripes in R. Phase-space 
trajectories confirm this picture, whereby the chaotic behavior is signalled by the absence of sharp closed trajectories. 


sate with a finite “phase-rigidity”, i.e. a finite energy 
cost oc f dx(d x (p ) 2 of creating a phase gradient in the 
wave function 'ip = \^ \ exp (iip), which is associated with 
a finite flow velocity of the superfluid v = hd x (p/m [49] . 

As illustrated in Fig. [4j the driven BEC coupled to 
the lossy cavity, once entering the nonlinear oscillating 
regime, shows large phase gradients and thus accumu¬ 
lates the corresponding energy. The most efficient way 
for the superfluid to get rid of this extra energy is to con¬ 
vert it into nonlinear dispersive waves through the pro¬ 
cess of phase-slippage ®t In one dimension these waves 
are solitons, for which phase-slips take place as described 
in Fig. [4] C). We observe indeed phase-slippage as the 
system enters the dynamically unstable regime. During 
the limit-cycles, phase slips take place periodically and 
are synchronised with the oscillations. The rate of slip¬ 
page is slow (see Fig. [4] D)) and the solitons are nucle¬ 
ated always at the same position, shown in Fig. @A1). 


On the contrary, once entered the chaotic regime, phase- 
slippage takes place at an orders of magnitude faster rate 
and the nucleation of solitons becomes irregular in space 
and time. 

Interestingly, this fast proliferation of solitons does not 
eventually destroy the BEC as it usually happens in an 
isolated superfluid. In the latter case the large num¬ 
ber of solitons causes the phase of the wavefunction to 
change fast and over short lenghscales, thereby destroy¬ 
ing phase coherence. In addition, solitons contain a large 
kinetic energy which can be converted into the thermal 
component. However, in our driven/dissipative system 
the proliferation of solitons is counteracted by light scat¬ 
tering processes transferring energy into the cavity field 
and expelling it through photon losses. Indeed, we oberve 
that even though solitons are produced at a fast rate, the 
BEC kinetic energy = f dxh 2 \d x ip\ 2 /2m, after an ini¬ 
tial increase, fluctuates around a finite value of the order 
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FIG. 4. Phase-slippage through soliton nucleation for a finite value of the atom-atom interaction g aa — l.OhuJn- Panels A) 
and B) show the condensate density (color-scale) as a function of position and time (A1 and Bl) as well as the time-evolution 
of the cavity amplitude (A2 and B2) for 77 = 12.0o;r, which corresponds to a limit-cycle dynamics (A) and 77 = 22.0o;r, 
which corresponds to chaotic dynamics (B). In A),B), white circles indicate the space-time coordinate of phase-slips, i.e. phase 
singularities in the form of solitons having a zero-density core together with a 7 r phase difference across it. The dynamics of a 
single phase-slip is illustrated in panels Cl) to C3) for 77 = 22.0cc;r. The phase singularity is present in C2) as a 7r-phase jump 
localized at the dark (n(x) = 0) soliton position. The latter has a size of the order of the healing length £ = l/y/2mg c n ~ 0.1A. 
The phase-slip consists in changing the sign of the phase jump n —)> —tv (they are equivalent at the point n{x) — 0). Subsequently, 
in C3) the soliton increases its minimum density from zero and moves away, carrying the energy subtracted from the initial 
phase gradient in Cl). Panel D) shows the total number of phase-slips N s {t) (solid line) generated up time t together with the 
condensate kinetic energy E k (dotted line) as a function of time, for 77 = 12.CL;r (in blue) and 77 = 22.0o;r (in red). While in the 
limit-cycle regime phase-slips are generated following the slow periodic oscillations (two phase-slips every second period in the 
case shown), the rate of slippage in the chaotic regime is orders of magnitude faster. Still, soliton proliferation is counteracted 
by dissipation through cavity losses so that in the chaotic regime the kinetic energy intially grows but then flucuates around a 
constant value. 


of the recoil energy /icjr, as shown in Fig. [4j 


V. CONCLUSIONS AND OUTLOOK 

Surprisingly, dynamical self-ordering of a BEC inside 
a Fabry-Perot cavity appears also in the blue detuned 
regime but with decisive new properties. While collec¬ 
tive scattering is enhanced by ordering it is inhibited by 
repulsion from field antinodes. This competing effects 
generate self-sustained large amplitude nonlinear oscilla¬ 
tions, which reconcile ordering and superradiant scatter¬ 
ing averaged over an oscillation cycle. For stronger pump 
the oscillations exhibit period doubling and evolve into 
chaotic behavior. This complex dynamics takes place 
without any external modulation or feedback and is an 
intrinsic property of the driven dissipative nature of the 


dynamics. The nonlinear atom-field dynamics here is ac¬ 
companied by phase-slippage, which comes about due to 
the superfluid nature of the interacting BEC. In one di¬ 
mension phase-slippage nucleates solitons whose prolif¬ 
eration in the chaotic regime is counteracted by cavity 
cooling. 

In two spatial dimensions vortex nucleation is instead 
expected. It would be interesting to study how their 
proliferation and dynamics is affected by light scattering 
into the lossy cavity. In this spirit, the extension of the 
present study to d = 2 should open up new directions in 
the field of quantum turbulence m in driven/dissipative 
condensates [51 . 

We finally point out that the system studied here is al¬ 
ready available experimentally since it involves a change 
from red- to blue-detuning of the driving laser with re¬ 
spect to the atomic resonance in the setups used in the 
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